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The original idea to show the spacetime geometry using few geodesies was devel- 
oped by Johnson and Ruffini (1974). We used this idea to interpret the observa- 
tional data for rotating BH's. We developed the imitation approach to simulate 
a propagation of radiation near BH's. An important problem for this approach is 
the diagnostics of a blacli hole metric using X-ray observational data of the iron 
Ka-line. Observations of Seyfert galaxies in X-ray region reveal the broad emissi- 
ion lines in their spectra, which can arise in inner parts of accretion disks, where 
the effects of General Relativity (GR) must be counted. A spectrum of a solitary 
emission line (the i^^-line of iron, for example) of a hot spot in Kerr accretion 
disk is simulated, depending on the radial coordinate r and the angular momen- 
tum a = J/M of a black hole, under the assumption of an equatorial circular 
motion of a hot spot. Using results of numerical simulations it is shown that the 
characteristic two-peak line profile with the sharp edges arises at a large distance, 
(about r f5! (3 — 10)rg). The inner regions emit the line, which is observed with 
one maximum and extremely broad red wing. High accuracy future spectral ob- 
servations, being carried out, could detect the angular momentum a of the black 
hole. We analyzed the different parameters of problems on the observable shape of 
this line and discussed some possible kinds of these shapes. The total number of 
geodesies is about 10^ (to simulate possible shapes of the ii'c-line), so the number 
is great enough, especially in comparison with few geodesies in the original paper 
by Johnson and Ruffini (1974). 

It is clear that an analysis of geodesies gives the direct way to investigate a met- 
ric. In particular, Johnson and Ruffinff calculated geodesies near a Kerr black hole. 
These geodesies were very beautiful and they were used to create the known ICRA 
emblem. However, there is a problem to recognize such shapes of these geodesies 
near black holes, because practically we use something like photon geodesies to 
look at the geodesies near black holes, but the photon geodesies are bent by a very 
complicated way, so it is extremely difhcult to reconstruct the shapes of geodesies 
near black holes. Below we will discuss how to extract an information about a 
metric from an analysis of geodesies around black holes. 

The general status of bl ack ho les described in a number of papers (see, for 
example the following 

paperj™ 

and references therein). As it was emphasized 
in these reviews the most solid evidence for an existence of black holes comes 
from observations of some Seyfert galaxies because we need a strong gravitational 
field approximation to interpret these observational data, so probably we observe 
manifestations radiation processes from the vicinity of the black hole horizon (these 
regions are located inside the Schwarzschild black hole horizon, but outside the 
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Kerr black hole horizon, thus we should conclude that we have manifestations of 
rotational black holes). 

Recent observations of Seyfert galaxies in X-ray band I I ^ 1 reveal the 

existence of wide iron line (6.4 keV) in their spectra along with a number 
of o ther weaker lin es (No X, Si XIII,XIV, S XIV-XVI, Ar XVII,XVIII, Ca XIX, 
etc) III ( | o | ij | iU | i i | rpj^g jjj-^g width corresponds to the velocity of the matter mo- 
tion of tens of thousands kilometers per second, reaching the maximum value 
V « 80000 " 100000 km/s for the galaxy MCG-6-30-1^ and v « 48000 km/s 
for MCG-5-23-16.C3 In some cases the line has characteristic two-peak profile with 
a high "blue" maximum and the lo w "red " one and the long red wing, which grad- 
ually drops to the background level.l^^i^ 

To simulate these shapes of the spectral lines we choose a minimal number of 
assumptions. We used the numerical approach based on the method, described 

Many astrophysical processes, where the great energy release is observed, are 
assumed to be connected with the black holes. Because the main part of the astro- 
nomical objects, such as the stars and galaxies, possesses the proper rotation, then 
there are no doubts that the black holes, both stellar and supermassive, possess the 
intrinsic proper rotation too. Therefore we consider an emission of monoenergetic 
quanta near a Kerr black hole. 

The large amount of observational data requires its comprehension, theoreti- 
cal simulation and interpretation. The numerical simulations of t he a ccretion disk 
spectrum under GR assumptions has been reported in the paper.EOI The observa- 
tional manifestations of GR effects are considered in X-ray binarieJ^. Different 
physical models of the origin of a br oad ern ission iron Ka line in the nuclei of Seyfert 
galaxies are analyzed in the papers.^ES 

The stationary black holes are described by the Kerr metricPSI 

ds^ = -4 (dt - asin^ ed<pf + ^ [(r2 + a") d<f> - adtf + ^dr^ + p'dO^ (1) 



where 

p2 = + cos^ e, (2) 

A = - 2Mr + c? . (3) 

The photons trajectories can be described by the standard equations of 
geodesies: 

Sx^ , dx^ dx'- , , 

where are the Christoffel symbols. The equations geodesies however can be 
simplified if we will use the complete set of the first integrals which were found by 
Carter*^ E = pt is the particle energy at infinity, = is z-component of its 
angular momentum, m = pip^ is the particle mass and Q is the Carter's separation 
constant: 

Q=pI + cos^ [a^ (m^ - E^) + Ll/sin^ 0] . (5) 
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As shown by ZakharovpEl the equations of photon motion can be reduced to 

+ + (6) 

(7) 

= 2r3 + {a^ -e-v)r+{a-0+ V, (8) 

(9) 

a2sin6'^, (10) 
r^ + a^-^a), (11) 



dt' 




da 


= —a 


dr 




da 


= ri, 


dri 
da 


= 2r' 


dO 




da 




dOi 




da 


= cos 


d<j) 
da 





sin^6' 



sin^ ej A 

where rj = Q/KPE^ and ^ = Lz/ME are the Chandrasekhar's const ant J^', which 
should be derived from the initial conditions in the disk plane; r and a are the 
appropriate dimensionless radial variable and constant of rotation respectively. The 
system (j^- (|ll|l has also two integrals, 

ti=rl-r^ ~ (a^ - _ ^2 _ 2 \[a - £,f + r/1 r + a^r] = 0, (12) 
£2 = 6*? - 77 - cos^ e (a^ - = 0, (13) 



which can be used for the precision control. This method differs from the approach 
which was developed in papers. "^^ ■^^ 

First, one should define the Chandrasekhar's constants for each quantum and 
then integrate the system (|5|)- (|llll to either the infinity or the events horizon, de- 
pending on the constants values. 

We assume that the hot ring emits quanta which are distributed by isotropic 
way in its local frame. The simulations are based on the trajectory classification, 
depending on the Chandrasekhar's constants ^-—'^^ ' 

The simulated spectrum of a hot ring for a = 0.9, 9 ~ 60° and different radius 
values is shown in Fig. ^ The proper quantum energy (in co-moving frame) is 
set to unity. The observer at infinity registers the characteristic two-peak profile, 
where the "blue" peak is higher than the "red" one and the center is shifted to the 
left. Some spectrum "oscillations" near its minimum is explained by pure statistical 
reasons and has no the physical nature. 

As far as the radius diminishes the spectrum is enhanced, i.e. increases the 
residual between the maximum and minimum quanta energy, registered by a distant 
observer. For example, for a = 0.9, r = 1.2 Vg and 9 = 60°, where rg has its 
standard form Tg = 2kM/(? , i.e. in the vicinity of the marginally stable orbit, the 
quanta, flown out to the distant observer, may differ 5 times in their energy. The 
red maximum decreases its height with diminishing the radius and at r < 2 is 
almost undistinguished. It is interesting to note that the spectrum has very sharp 
edges, both red and blue. Thus, for a = 0.9, r = 3 rg, 9 = 60° the distant observer 
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Figure 1. Spectrum of the hot ring for a = 0.9, 9 = 60° and different radial coordinates. The 
marginally stable orbit lays at r = 1.16 rg. 



has registered 1433 quanta of 20417 emitted by isotropic way; 127 of them (w 
9%) drop to the interval 1.184 < E < 1.202 (blue maximum) and 43 quanta drop 
to 0.525 < E < 0.533 (red maximum), whereas no one quantum has the energy 
E < 0.518 or £; > 1.236. 

A spectrum of a hot ring for a — 0.9, r — 1.5 rg and different 9 values is shown 
on Fig. 121 The spectrum for 9 — 60° and the same a and r values is included 
in Fig. n and should be added to the current figure too. As it follows from the 
figure, the spectrum critically depends on the disk inclination angle. For large 9 
values, when the line of sight slips almost along the disk plane, the spectrum is 
strongly stretched, its red maximum is essentially absent, but the blue one appears 
narrow and very high. The red wing is strongly stretched because of the Doppler 
effect, so that the observer registers the quanta with 5 times energy difference. As 
far as the 9 angle diminishes the spectrum grows narrow and changes the shape: 
its red maximum first appears and then gradually increases its height. At 9 = 0° 
both maxima merge to each other and the spectrum looks like the (5-function. It 
is evident because all the points of the emitting ring are equal in their conditions 
with respect to the observer. The frequency of registered quanta in that case is 2 
time lower than the frequency of the emitted ones. A fall in frequency consists here 
in two effects, acting in the same direction: the transversal Doppler effect and the 
gravitational red shift. 

The strong variability of Seyfert galaxies in X-ray does not contradict the as- 
sumption, that we observe the emission of the hot rings from the inner region of 
accretion disk, which can decay or grow dim, going towards a horizon as time 



ZAKHREPI: submitted to World Scientific on February 2, 2008 4 



rl 1.0 n 



:0.6 =0.6 



:0.2 0.2 




Figure 2. Spectrum of a hot ring for a = 0.9, r = 1.5 Tg and different 6 angles. 



passes. The spectrum dynamics is understood qualitatively by reference to Fig. ^ 
considered sequentially from top to bottom. 

To analyze an influence of a disk width on the shapes of the line we consider the 
case of a wide accretion disk a nd it w as shown that the shape of the spectral line 
retains its type with two 

peakJI3S2 (see FigOJ. It is noted that the inner parts 
give the essential contribution into red wing of spectrum. 

It is known that the standard disk models (like, for example, Shakura - Sun- 
yaev and Novikov - Thorne disk models) hardly ever could be used to describe 
temperature distributions in accretion disks of Seyfert galaxies, however to show 
an influence of a temperature distribution on the spectral line shapes we use the 
standard disk model as a template. Fig. ^ demonstrates the shape of emitted 
monochromatic line in Sc hwar zschild black hole fleld with temperature distributed 
according to a-disk model.l^ 

De tails of c omputations and a full list of references could be found in 
papers.^^l^ An application of such approach to estimate magnetic fields in 
AGNs and microquasars is described in 

details.EHl 
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Figure 3. The spectral line shapes for different 6 angles. The emitting region is the wide ring and 
its inner boundary is the last stable orbit (for rotational parameter a = 0.9 this r-value is equal 
to r = l.Wvg), its outer boundary corresponds to r = lOr^. 
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Figure 4. The spectral line shapes corresponding to an accretion disk with outer and inner radii 
{fin = 3rg and rout = 10 rg) in Schwarzschild black hole field for different position angles of a 
distant observer. The temperature is distributed a<;cording to the ShaJfura - Sunyaev model. 
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